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Abstract. A new class of cost functionals for optimal control of quantum systems which 
produces controls which are sparse in frequency and smooth in time is proposed. This is 
achieved by penalizing a suitable time-frequency representation of the control field, rather 
than the control field itself, and by employing norms which are of or measure form 
with respect to frequency but smooth with respect to time. 

We prove existence of optimal controls for the resulting nonsmooth optimization prob¬ 
lem, derive necessary optimality conditions, and rigorously establish the frequency-sparsity 
of the optimizers. More precisely, we show that the time-frequency representation of the 
control field, which a priori admits a continuum of frequencies, is supported on only 
finitely many frequencies. These results cover important systems of physical interest, in¬ 
cluding (infinite-dimensional) Schrodinger dynamics on multiple potential energy surfaces 
as arising in laser control of chemical reactions. Numerical simulations confirm that the 
optimal controls, unlike those obtained with the usual Lfi costs, concentrate on just a few 
frequencies, even in the infinite-dimensional case of laser-controlled chemical reactions. 


1 Introduction 

This paper is motivated by application problems of current interest in quantum con¬ 
trol nn, which range from steering chemical reactions [2] over creating excited or ionized 
states |12j to faithfully storing and manipulating bits of quantum information m- 

We propose a new class of cost functionals for the optimal control of quantum systems 
which result in controls with a sparse time-frequency structure. This is achieved via two 
key ideas. 

First, we do not penalize the time profile of the field amplitude but a suitable time- 
frequency representation of it. While such representations are a familiar tool to interpreting 
or analyzing a given field in control and signal analysis, they here acquire center stage 
already in the design of the controls. 

Second, we build upon recent advances in the optimal control theory of elliptic and 
parabolic systems related to the basic idea [39l EH El El HH of sparsity-enhancing or 
measure-norm costs. More specifically, we build upon the idea of function-valued measures 
to achieve directional sparsity in parabolic control m- The novelty as compared to the 
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latter advances is two-fold: first, in quantum control, unlike in parabolic control, the target 
for sparsity should not be the field amplitude, but its frequency structure, and second, 
one is dealing with a bilinear instead of a linear control problem. 

These ideas result in constrained non-smooth optimization problems of the form 

Minimize -|-Q;||ri||_A 4 over controls tt: x [0, T] —)• C (1) 

subject to 

idt'ip = {Ho + {Bu){t)Hi)'ip, •!/)(0) = ■i/’o- (2) 

Here Hq is the Hamiltonian of the quantum system, H is a ‘synthesis operator’ which 
assembles the control held from a time-frequency representation u{uj,t), and || • \\m is an 
or measure norm with respect to frequency but a smoothness-promoting norm with 
respect to time. A prototypical choice is the time-frequency synthesis operator 

{Bu){t) = [ n(a;, t) dw, (3) 

Jn 

where H C M is a region of admissible frequencies, and 


Mm = / lk(w,-)lbi(o,T)da;. 

Jn 


( 4 ) 


Note that the control, a priori, can use the whole available continuum of frequencies, with 
each frequency possessing its own time prohle. A main result of this paper (Theorem |4.8[ ) 
is that the optimizers utilize only finitely many frequencies, even when the quantum 
dynamics is a full inhnite-dimensional Schrodinger dynamics on multiple potential energy 
surfaces (Ex. 2.3 and Section 5.3). Equations (|^-Q replace the standard approach in 
quantum optimal control initiated by |26j to penalize just the Lfi or norm of the held 
amplitude (see [IZIETIIIS] for mathematical results). 

Numerical simulations presented in Section below illustrate that the optimizers con¬ 
centrate on just a few frequencies, even when the quantum dynamics is a full inhnite- 
dimensional Schrodinger dynamics. Thus our controls share an important feature of laser 
pulses designed by experimentalists. 

The plan of this paper is as follows. In the next section we introduce three important 
physical examples of controlled quantum dynamical systems which motivated this work 
and to which our sparsity results apply. These examples also serve to recall basic features 
of the coupling operators such as ‘forbidden transitions’ and the oscillatory nature of 
quantum controls. In Section]^ we introduce our measure-norm sparsity-enhancing costs 
within a general functional-analytic framework, and give several examples. In particular, 
the choices 0-S lead to frequency-sparsity with global time prohles, and appropriate 
modihcations lead to frequency-sparsity with local time prohles. Section is devoted to 
the mathematical analysis of the non-smooth optimal control problem We recall 

the relevant well known results on the existence of dynamics, establish existence of opti¬ 
mal controls, derive necessary optimality conditions, and prove that optimal controls are 
supported on only hnitely many frequencies. Finally, in Section we numerically calcu¬ 
late optimal controls and compare them to those obtained from the usual Lfi penalization 
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of the field amplitude. Specifically, we present a 3-level example which arises in atomic 
excitation problems, and an example of Schrodinger dynamics on two potential energy 
surfaces as arising in laser-controlled chemical reaction dynamics. 


2 Quantum dynamics with controls: physical examples 

The evolution equations in quantum control problems typically have the structure 

V'(t), (5) 

where the state '0(t) belongs to some Hilbert space H, Hq and the Hi are (bounded or 
unbounded) self-adjoint operators on T-l, and the V£{t) are real-valued scalar amplitudes 
of components of applied electric or magnetic fields. The operator Hq is the Hamiltonian 
of the system in the absence of fields, and the operators describe the system-held 
coupling. 

A general mathematical reference for hnite-dimensional problems of form ([^ is [10]. In 
the mathematical control theory literature, general aspects of inhnite-dimensional prob¬ 
lems of the above form such as existence of optimal controls with or penalization 
of the held have been previously treated HZIEE]. Here our goal is to mathematically 
understand the highly oscillatory nature of optimal controls, unfamiliar from elliptic and 
parabolic control problems but arising in the physics of Bohr frequencies, and to develop 
novel penalty terms to simplify this oscillatory structure. 

We note that equation ([^ already contains two important approximations which are 
valid in many situations of interest. First, quantum huctuations of the held amplitudes 
can be neglected, that is to say we are dealing with classical helds and do not need to move 
to the much more complicated framework of quantum held theory. Second, the spatial 
wavelength of the applied helds is much larger than the localization length of the state 
so that it is sufficient to assume that the held strengths depend on time only. This 
is often called ‘dipole approximation’, see e.g. m- 

We now give three examples for Q of physical interest. The hrst one has been included 
to recall the quantum mechanical meaning of oscillatory controls which forms the starting 
point for the time-frequency approach developed here; the other two will be used in the 
simulations in this paper. 

Example 2.1. (Spin of a spin 1/2 particle in a magnetic held) This is the simplest control 
system of physical interest. It arises as a basic example in NMR, and more recently as 
a model of a single qubit in quantum information theory (see [3l| for a recent careful 
experimental realization of this system and |5| for rigorous mathematical results). It 
already exhibits surprisingly many features of complex systems. The spin at time t is a 
unit vector in the Hilbert space H = C'^. The general evolution equation of a spin in a 
time-dependent magnetic held H: M —>■ is 

3 

idtifit) = ■ Si){t) = -7 ^ Ba{t)Sa ifit) 

e=i 
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where the component operators Sa of the spin operator S are given by h /2 times the Pauli 
matrices, i.e. in atomic units {h = 1 ) 


Si 




S2 



S3 



The factor 7 depends on the type of particle (electron, proton, neutron, nucleus) and can 
be positive or negative. A typical control problem in NMR consists of taking B 3 time- 
independent and comparatively large, and Bi and B 2 as time-dependent controls which 
are small. This is a system of form ([^, with Hq = —’-yB^S^. Denoting the two eigenvalues 
of Hq by El, E 2 , this system can be written in the compact form 


idtii 









with complex-valued control v{t) = —2'y[Bi{t) +iB 2 {t)), and v* its complex conjugate. 
The basic case of a time-harmonic control field v{t) = is exactly soluble [27j. This 

allows one to understand mathematically the emergence of oscillatory controls and Bohr 
frequencies. The Bohr frequency of a transition between two quantum states is the eigen¬ 
value difference. The time-harmonic control with this frequency, when applied over a time 
window of suitable length, achieves a 100 % transfer; and it is the only time-harmonic 
control which achieves a 100 % transfer m- For more general quantum systems relations 
of this kind typically only hold asymptotically [7]. 


Example 2.2. (Electronic states of atoms in laser fields) A standard reference in the 
physics literature is |33j . Consider an atom with N electrons of charge e = — 1 and a nucleus 
of charge Z = N clamped at the origin. The electronic state of the atom is described by 
a function belonging to the Hilbert space = {'!/> G L^((M^ x Z 2 ) ) | '(/> antisymmetric}. 
That is to say electronic states are functions = ^^{xi, si,... ,xn, sn) which depend on 
the position coordinates Xi G and the spin coordinates s* G Z 2 of all the electrons. An 
applied electric field can be described by a function E: M —)> with E{t) denoting the 

electric field vector at time t. (Here quantum fluctuations of the field as well as its spatial 
dependence are neglected, as discussed above.) The overall evolution equation is 


idt'ip{t) 


--A + V{xi,... ,xn) - E{t) ■ D{xi,... ,xn) ) -ipit), 


( 6 ) 


where A is the Laplacian on , the many-body Coulomb potential V is given by 
V{xi,... ,xn) = + Yli<i<j<N^/\^i ~ ^ is the dipole operator 

D{xi ,..., xn) = exi- This has the form ([^, as is immediate by denoting the com¬ 
ponents of E{t) and Xi with respect to some orthonormal basis of by Ei{t) and xa 
(.^ = 1, 2, 3), and letting Hi = Because the high-dimensional Schrodinger equa¬ 

tion Q cannot be simulated in practice, the infinite-dimensional state equation ([^ is 
often replaced by projecting onto finitely many eigenstates V'l, • • • j of Hq = — ^A -|- V, 
i.e. ip{t) Ri ai(t)V’i -!-••• + ad{t)'ijjd and neglecting the coupling with the rest of the system. 
Assuming for simplicity that the field is unidirectional, i.e. E(t) = v{t)EQ for some unit 
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vector Eq and a scalar amplitude v, and denoting the eigenvalues of Hq corresponding to 
the states tpi,... by Ei,..., E^, this yields 
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with the coupling matrix elements 

fJ'mn — (V’mj Eq • D Ipn) ■ 


( 7 ) 


( 8 ) 


We remark that the off-diagonal element /i^n vanish whenever ipm and 'ijjn have the same 
parity. This is a simple example of a forbidden transition. 


Example 2.3. (Laser-guided chemical reactions) Our last model is of central interest in 
photochemistry, but to our knowledge has not hitherto been considered at all in the math¬ 
ematical literature. Consider a molecule with M nuclei at positions R = (i2i,..., Rm) £ 
and N electrons at positions xi,...,XAr with spins si,...,SAr. The state of the 
molecule at time t is described by a wave function T G x (M^ x Z 2 ) ), i.e. 

T = '^{R, xi, si,..., X]sf, si\f). The evolution equation is given by 




/ M 

E 

v.a=l 


Ar, + + E{t) • 


2 mn 




with electronic Hamiltonian = — b/S.+V^^\xi, ■ ■ ■ ■, xn), potential ... ,X]\f) = 

~ Sa=l ^al\xi ~ Ra\ + Yll<i<j<N 7/1^* ~ Xj\ + Yll<a<^<M ~ and 

dipole operator D^^^{xi,... ,X]\[) = ~ 'l2a=i^^oiRa- Here, e is the electronic 

charge (—1 in atomic units), —eZ^ are the nuclear charges {+Za in atomic units), and ma 
are the nuclear masses. A careful mathematical account in the absence of control helds 
and for smooth interaction potentials is given in [35]. 

The typical situation in laser control of chemical reactions is the following, see Figure 
The system starts in a stationary state of electrons and nuclei. The laser then induces a 
transition to a different electronic state. As a result the nuclei now see a different potential 
energy surface with respect to which they are no longer in equilibrium; for instance the 
new surface may no longer contain a barrier to a desired target position. Once the nuclei 
have moved barrier-free to the target position, the laser induces a transition back to the 
original surface so as to also put the electrons in the target state. Mathematically, this 
situation can be modelled by generalizing the ansatz for electronic state in Example |2.2| 
to an ansatz of the joint wave function of electrons and nuclei. Confining ourselves for 
simplicity to two electronic states, one assumes 


'^{R,Xi,Si,...,XN,SN,t) 


^i{R,t)il}\ \xi,si,...,xn,sn) 


-h $ 2 ( 7 ?, t )'02 {xi,si,...,xn-,sn), 
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Figure 1: Schematic representation of laser-controlled chemical reaction dynamics. The 
nuclei of a molecule move on different potential energy surfaces depending on 
the electronic state, and the laser induces transitions between these states. Blue: 
Potential energy surfaces. Magenta: Initial wave function of the nuclei. Cyan: 
Target region. 


where <I>i, <^2 G are nuclear wave functions, ||‘hi|P + ||‘h 2 |P = I, and V’ 2 ^^ 

/ ry\ 

are normalized eigenstates of the electronic Hamiltonian . This leads to the following 
Schrodinger equation in the Hilbert space Ti = C^) 



+ yE{t) • fii2{Rr m • fi22{R)j) [^ 2 )^ ^ ’ 

with the dipole moment functions 

= ( 10 ) 


3 Cost functionals and functional analytic setting 

To identify control helds vi: [0, T] —?• M which achieve a suitable goal, such as transfer of 
the system from one eigenstate of Hq to another, and which are convenient to implement 
experimentally, one typically minimizes a cost functional which promotes goal achievement 
and penalizes unsuitable controls m- The theoretical work up to now has focused on 
cost functionals like the norm of the control field (or variants thereof like the or 
norm). This approach has been very successful in achieving a wide spectrum of control 
goals. But it has been remarked that “[the] resulting optimized electric fields are often 
very complex, thus it is nearly impossible to understand the underlying processes involved 
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in the observed control” |16] . We may therefore ask the question: Why (or variants 
thereof)? 

We will argue that using certain more sophisticated costs leads to “simpler” controls. 
The control helds are sparse in frequency, picking out the system’s Bohr frequencies with¬ 
out any special ansatz, while at the same time having slowly varying amplitude envelopes, 
thereby sharing an important feature of laser pulses designed by experimentalists which 
are commonly and successfully used in the laboratory. 

We first state the general class of costs we propose, including its functional analytic 
setting. Subsequently we give examples. The simplest example of such a cost which 
yields frequency-sparse controls, and which motivated the general setting, is described in 
Example 3.1 below. 


3.1 General setting 

For the quantum system 


L 

idtifit) = [Ho + ^ Bu{t)iHi)'tjj{t), V’(O) = V’o (H) 

i=i 

we consider the optimal control problem 

min J{'ip,u) subject to pTj) (12) 

ijj, U '-' 

where 

= ]^{'ii){T),0'ip{T))'n + a\\u\\M{a-M)- ( 13 ) 

The functional J consists of the term ^{i/;{T),0'ip{T))'^ which describes the expectation 
value that needs to be minimized, and a cost term of the form a||rr||x(n;W) which contains 
a measure norm explained below and forces the solution to be sparse in a suitable sense. 

We make the following very general functional-analytic assumptions on the operators 
and fields appearing in (|ll[)-([l2)) . These assumptions cover all the physical examples from 
Section [2l 

1. Dynamics. Assume that the Hamiltonian Hq is a self-adjoint operator on a Hilbert 

space Ti with domain T>{Hq). The initial condition ■i/’o may be any element of %. The 
coupling operators are assumed to be bounded self-adjoint operators on T-L. We use the 
vector operator notation v H = vi{-iHi), v H* = Ya^i vi{-iHi)*, {xi,HX2)h = 

((Xi, -iHiX2)'H)i=i e for H = {-iHi,..., -iHi), v and xi,X 2 G %. 

For the admissible class of controls u and control operators Bi specified below, we will 
show that eq. (0 possesses a unique mild solution in the state space of continuous paths 
in the Hilbert space, C'([0, T]; H). 

2. Target constraint. The observable O specifying the target constraint can be 
any bounded self-adjoint operator on B. Typically O is the orthogonal projection onto 
the subspace we want to reach. If O is a projection, the target constraint contribution 
^{'if{T),Oif{T))'^ to the cost lies in the interval [0,1/2]. The value 0 corresponds to a 
100% achievement of the control objective, and the value 0.5 to a 0% achievement. 
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3. Control space, cost, measure norm. The control field u is assumed to belong 
to a measure space of form 

(14) 

where is a locally compact space (typically a closed interval of admissible frequencies), 
and U is a separable Hilbert space of time-dependent functions (admitting general sepa¬ 
rable Hilbert spaces can be useful to obtain nice optimality conditions, see the discussion 


after Example 3.2). The space (14) is the space of Borel measures u on H with values in 

|tt||_yK. The mass norm of the measure u is the second term in the 


U of finite mass norm 


cost functional J in (13). 


The space (14) is the dual of the separable space Cq{VL]U) of continuous functions on H 


with values in lA which can be uniformly approximated by functions with compact support. 
The duality pairing is given by 


{u,^)m,Co = / {u'{u}),ip{u;))ud\u\{uj) (15) 

Jn 

(see [24] i where u' is the Radon-Nikodym derivative of u with respect to the total variation 
measure |u| (see I2Q1 VII Thm. 4.1]). Note that the inner product in the integral is the 
Hilbert space inner product in U. This duality will be very useful. 

4. Control operator. The control operator is assumed to be a bounded linear 
operator 

B: JiA{Q]U) —)• T^( 0 , T; M^) for some 1 < p < oo. 

Moreover we assume that B has a bounded linear predual operator B*: L'^(0,T;M^) —?■ 
CoiQ'M), by which we mean a bounded linear operator such that (B* f, u)cn M = (f, Bu) rq rp 
tor ajl/einO,r;Mqa„dall.,eA 1 (ii,M)wLe 1 + 1 = 1 . 

Existence of a bounded linear predual operators implies the weak-=t:-weak(->i:) continuity 
of B. That is, weak-* convergence of Un to u in A4(H;^) implies weak convergence of Bun 
to Bu in LP{0, T ; M^) if 1 < p < oo, and weak-* convergence if p = oo. The case p = 1 has 
to be excluded in our framework since we will make use of the reflexivity of LP( 0 ,T;M^). 

Note that the operator B* depends on the Hilbert space structure of lA, see examples 
below. Since B* appears in the optimality system, the freedom to choose lA can be used 
to generate nice optimal controls. 

All spaces — possibly containing complex valued objects — are equipped with a real 
Banach or Hilbert space structure. For complex spaces this amounts to always using the 
real part of the scalar product, i.e. we take {ip,'ip)p to mean Re((p,particular, 
linear always means M-linear and the scalar product is real-valued and M-bilinear. 

3.2 Examples 

We now list some examples for choosing the frequency domain H, the Hilbert space lA of 
time-dependent functions, and the control operator B. The first example is the proto¬ 
type for generating sparse controls. It motivated the general functional analytic setting 
proposed above, and naturally incorporates physically relevant controls containing finitely 
many pulses of particular frequencies [il[ilEgE3En]. 






Example 3.1. (Two-scale synthesis, smooth functions of time) Here u will be a two- 
scale time-frequency representation of the laser field amplitude and B will generate the 
corresponding field. The control u can a priori contain a continuum of active frequencies, 
each with its own smooth envelope. This can be modeled mathematically as follows. Let 
n be a closed subset of M"*", U = Hq{0, T; C), and p = oo. For u G ILq (0, T; C)) we 

define B to be the two-scale synthesis operator 


(i?u)(t)=Re f duj. (16) 

Jn 

By approximation, the expression can be extended to measure-valued controls. This ex¬ 
tension is important in practice, because it allows sharp concentration on a small number 
of frequencies, and has the following mathematically rigorous integral representation: 


{Bu){t) = Re [ u (i\u\{oj), 

Jn 


(17) 


where u' is the Radon-Nikodym derivative of u with respect to |tt|. 

We remark that this setting naturally contains the physically motivated finite-dimensional 
ansatz spaces of [D EH EH [m in which a finite number of distinct frequencies can be 
switched on or off by few-parameter modulation functions: the field v{t) = cos{ujjt+ 

4>j) corresponds to Bu with B as in (© and u{Ld,t) = S{uj — L 0 j)e'^’^^Vj{t). We allow 

U to contain complex-valued functions. This allows phase shifts in the different frequencies 
without leaving the linear setting. 

The predual operator B* is the solution operator of the second-order boundary value 
problem 


dt'^ 


u{uj, t) 


fit)e 


—iojt 


U 


t=0 


= U 


t=T 


= 0 , 


(18) 


i.e. B*f = u. The equations ( [T^ are not coupled for different w. The operator B* is 
continuous and has the additional regularity B* f £ <70(^1; n Hq). Here it is important 
that n is closed. Otherwise B might not be weak-=t=-weak(-*) continuous. 


Example 3.2. (Gabor synthesis, functions of time) In this example we design the 
control operator B and the Hilbert space so that the predual operator B* has a nice 
form. We take C M"*", U = L^(0, T; C), and define B by 


(Ru)(t) = Re f [ k{s, t)u {u!, s) ds d\u\{uj) , (19) 

Jn Jo 

where k: [0,T]^ —)• M is a smooth symmetric kernel. Roughly, this operator corresponds 
to a pre-processing of the envelopes, with only smoothed envelopes entering the equation. 
The predual operator B* becomes the Gabor transformation 


{B*f){u,t)= f /(s)A:(t,s)e 
Jo 

which is a time-frequency representation of the control field. 
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Example |3.l| above with W = Hq leads to a time-frequency representation u = B* f which 
is global in time. That is, it has a global window equal to the Green’s function G of the 


one-dimensional boundary value problem (18), hence (i?*/)(a;, t) = L f(s)G(t,s}e 


On the other hand. Example |3.2| leads to a time-frequency representation which is local 
in time, but the definition of B is somewhat complicated. Note that Example 3.2 


can 


be reformulated in terms of the control operator from Example |3.1| by using a different 
control space. To see this, let K : L^(0,T) —)• L?‘{Q,T) be the compact operator given by 
convolution with a Gaussian kernel k, {Kf){t) = k{t, s)f{s) ds. Then K is injective and 
self-adjoint and has an unbounded inverse A: —)• L^. Define Uk '■= B{A) with the 

induced scalar product {u,v)u = {Au, Av) 12 . Then Uk is a Hilbert space and the predual 
operator B* of B is (B*f)(u},t) = (s). This construction also works for 

window functions other than a Gaussian. The equivalence to the choice U = L^(0,r;C), 
with H as in (19) follows since the dual of the map X : C'o(D; L^(0, T; C)) — >■ Co{B]ly(k) 


defined by {X(p){u}) = Kip{uj) is an isometric isomorphism between the control spaces that 
preserves the image under the corresponding control operators. 

Our next example shows that our setting also covers interesting cases when U does not 
contain time-dependent functions. 

Example 3.3. (Fourier synthesis, constant functions of time) Let H C ZY = C, and 
let B be the Fourier synthesis operator, that is to say 


{Bu){t) = Re / u{uj)e 

Jn 


iu)t, 


(ia;( = Re / d|u|(a; 


The function u{oj) here can be viewed as a time-frequency representation u{u,t) which 
is independent of time t. The predual operator is, up to a constant factor, the Fourier 
transform of functions restricted to [0,T], 

{B*f){u) = f{u)= r/(t)e—M l 
Jo 

An alternative approach to achieving sparsity of a time-global frequency decomposition 
via an L? cost combined with iterative ‘frequency sifting’ is given in [29] . 

Example 3.4. (Time-frequency Gabor synthesis) Let H C M"*" x [0, T] be a subset of 
time-frequency space, L/ = C, and 


(Ru)(t) = Re [ u'{uj,s)guj,sit)d\u\{uj,t) 

Jn 

for the ansatz function gu}^s{t) = e 2 ^^ £qj. gQj^g a > 0. This defines a suitable 

extension of the control operator from Example |3.2| to measures in both time and frequency. 
With this control operator, each Dirac measure u = corresponds to a Gaussian wave 
packet centered at time t with frequency uj. The predual of the control operator is given 

by 

fi'^,t('S)/(s)ds. 

Jo 
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These examples by no means exhaust our framework, but are meant to give an idea 
of its flexibility. We also remark that frequency constraints could easily be included by 
restricting fl, compare | 21 ) . 


4 Theory of the optimal control problem 


In this section we will study the optimal control problem We hrst show well-posedness 
of the problem. In contrast to m the main difficulty does not lie in the low regularity of 
the control since we assume sufficient smoothing of the control operator. It rather lies in 
the bilinearity of the state equation together with the low regularity of the state. Subse¬ 
quently we derive necessary optimality conditions. We shall show that the support of the 
optimal measure can be influenced. Theorem 4.6 is the natural analog of Theorem 2.12 
19). Differences arise due to the bilinearity of the equation and the non-trivial con- 


m 

trol operator. We shall also note interesting relationships between the choices for B and 
U. Throughout this section we will stay in the setting of mild solutions. This suggests 
to develop a derivation of the necessary optimality conditions which only requires inte¬ 
gral manipulations and no further regularity discussion for the primal and dual state are 
necessary. 


4.1 Existence and compactness properties for the state equation 

Here we recall well-known properties of the state equation (0 needed to prove the exis¬ 
tence of solutions to the optimal control problem. Throughout this subsection, for a given 
control field v G L^{0,T;'K^) we consider mild solutions to the state equation 0 i.e. 
functions ^|J which satisfy V' G C{[0,T];T-l) and 

ip{t) = + f G{t — s)v{s) ■ H'i/j{s) ds. (20) 

Jo 


Here G is the unitary group generated by the skew-adjoint operator —iHq and H = 
{—iHi)i. Existence of mild solutions as well as their differentiability in the direction of the 
field follow from [S) Thm. 2.5], and a convenient expression for the derivative is given in |22l 
Chapter 2]. The existence of mild solutions and their differentiability in the direction of 
the held are classical results |3l Theorem 2.5]. 


Proposition 4.1. Let v G L^(0,T;M^). Then (20) possesses a unique solution V' £ 
C{[0,T]-,'H). Furthermore, the mapping F: L^( 0 , T-M^) —> C{[0,T]-,'H) defined by F{v) = 
Tp is differentiable. The derivative is given by Ffiv){ 6 v) = where 


il)'{f)= [ G{t, s)Sv{s) ■ HtP{s) ds 
Jo 

with the evolution operator Q defined through Q[t, s)fi{s) = 'il>{t) for t > s. 


( 21 ) 


We remark that the evolution operators Q{t,s) are unitary. This can easily be shown 
by working in the rotating frame, i.e. by studying the function 1 1 —)• G*{t)'p{t) [28], but it 
will not be needed in the following. 
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One downside of the mild solution framework is that the typical compactness arguments 
based on additional spatial regularity cannot be used. Fortunately there is a 

powerful replacement that works in a very general setting. We will use the following 
compactness result from [3l Theorem 3.6]. 


Proposition 4.2. Let {vn)n 
corresponding solutions 'ipn of (20) satisfy ip. 
solution corresponding to v. 


he a sequence in L^(0,r;M^) such that Vn v. Then the 
■ Ip in C{[0,T];'H), where ip is the mild 


Remark. The results of Proposition |4.1| and 4.2 hold for general (Po-semi groups. The 
unitarity of the evolution operator is not used for the theoretical results in this work. 


4.2 Existence of optimal controls 

We will now prove the existence of solutions to the optimal control problem 


Theorem 4.3. There exists a global solution {ip,u) G C{[0,T]-,'H) x of (12). 


Proof. The result follows from standard reasoning in the calculus of variations. Let 
{ipn,Un) be a minimizing sequence, 


lim J{ipn, Un) = inf J{ip, u) 


( 22 ) 


Since O is bounded and a > 0, the sequence (un) is bounded in J\A(Q;IL}. There¬ 
fore, and because J\A(Q;IL) = C{Tl;U)* is the dual of a separable Banach space, the 
sequence {un)n has a weak-* convergent subsequence still denoted by {un)n with limit 
u G A4{fl;U). The weak-*-weak(-*) continuity of B implies that Bun converges to Bu 
weakly in T^(0, T; M.^) for some p > 1 and thus also for p = 1. By Proposition 4.2 the cor¬ 
responding sequence of states {ipn)n satisfies ipn{T) —>• ip{T). Thus the first summand of 
J converges, {ipn{T),Oipn{T)) —)• {ip{T),Oip) . The second summand of J is weak-* lower 
semi-continuous as it is a norm in a dual space. Thus we obtain lim„ J{ipn, Un) > J{ip, u). 
Together with (22) this implies the claim. □ 


4.3 Necessary optimality conditions 


For theoretical and numerical purposes we will use the reduced form of (12), 

minj(n), j{u) = J{ip{u),u). 


(23) 


Here ip{u) denotes the solution of for the control u. 

The reduced cost functional can be split into two parts. A nonlinear differentiable part 


fiv) = {iP{T),OiP{T))n 

with V = Bu, and a nondifferentiable convex part 

g{u) = a\\u\\M- 


(24) 
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In the next lemma we will see that the derivative of the differentiable part / is given by 

f{v) = 

where we recall from Section that {ip, 4’)'H denotes the real part of the inner product, 
and (/9 is the mild solution of the dual equation 

idtifit) = {Ho + v{t) ■ iH)ip{t), ip{T) = Oil){T). 

Using the evolution operator this can be rewritten as 


^{t) = g{T,trom. 


(25) 


Under suitable assumptions, the representation (24) of /' can be derived using a Lagrange 
functional approach, see 


is]. Since in our setting a variational formulation is not 
readily available, we will give a short proof in the present setting of mild solutions. 

Lemma 4.4. Let v, 6 v G L^(0, T;]R^), and let ip, ip' and cp be the corresponding solutions 
of (20), (21) and (25), respectively. Then the mapping f: L^{0,T;M.^) —^ M defined by 

f{v) = l{fi{T),Ofi{T))n 

is continuously differentiable with derivative 


f{v){Sv)= / 6 v{t) ■ {ip{t),Hfi{t))ndt. 

Jo 

Proof. The continuous differentiability of the state and the product rule give continuous 
differentiability of / and 

f'{v){5v) = {Ofi{T),fi'{T))n. 


Using (21) and (25) we obtain 


{Ofi{T),fi\T))n = {Ofi{T), / g{T,t)6v{t)-Hmdt)n 

Jo 

= [ 6vit) • {g{T, tyofiiT), Hm)n dt 
Jo 

= [ ^v{t) ■ {‘f{t),H'tp{t))ndt. 


□ 

We can now derive the following optimality condition. It is partially inspired the opti¬ 
mality condition derived in m Thm. 2.11] for a linear parabolic control problem of form 
dtip — Afi = u. 
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Proposition 4.5 (Optimality conditions). Let u he a local minimizer of problem (23), 
and let G C([0,T];'?^) be the corresponding solutions of (20) and (25) for the control 
field Bu. Then 

a\\u\\M = ( 26 ) 

and 

\\B*{^,Hf;)H\\co<a- (27) 

Proof. Since we can split our functional into a sum of a nonconvex and a nonsmooth part, 
the result can be deduced from the general differential calculus of Clarke [8] . Because this 
calculus is somewhat intricate and due to the importance of the optimality system, we 
prefer to give a more elementary proof. 

Let M be a local minimizer of problem (23) and let if and be the corresponding 

(28) 


solutions of (20) and 


. We first show the variational inequality 
9 {u) - f'{Bu){Bu - Bu) < g{u). 


Since u is locally optimal, we have 

^{j{u + h{u-u)) -jiu)) > 0 

for u G Ai{0,]hl) and h G (0,1) sufficiently small. Using the decomposition j = f o B + g 
and convexity of g, this implies 

i {f{Bu + h{Bu - Bu)) - f{Bu)) + g{u) - g{u) > 0. 


Since / is differentiable, taking the limit /i —)• 0 yields (28). 
Testing (28) with u = 0 and u = 2u gives 

g{u) + f'{Bu){Bu) = 0. 


Substituting (29) into (28) gives 


- f'{Bu){Bu) < g{u) 


(29) 

(30) 


for all u G M.{Vl]U). 


Using Lemma 4.4 on the derivative of /, equation 


gives 


9{u) = -{{(p,H'il;)n,Bu)L<i,Lp = -{B*{(p,H'il;)n,u)co,M 


which proves (26). From (30) we obtain 


-{B*{^p,H'tp)n,u)co,M < a\\u\\M- 

Testing this inequality with u = —6i^{B*{ip, Hfj)-}i){Lo) for some cj G yields 
\\iB*{p,Hf;)n)iu;)\\l < a\\iB*{p, Hf;)n)iu;)\\u 


which gives (27). 


□ 
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Remark. Proposition 4.5 provides only a necessary condition for local optimality. Due to 
the nonlinear structure of the problem (12), we expect that there also exist non-optimal 
critical points of j, as well as local optima that are not global. 

Proposition implies the following interesting restrictions on the support and direction 
of the optimal measure. 

Proposition 4.6 (Characterization of frequency support). Let u, and (p he as in Propo- 
sitionU.^ Then we have 


supp|u| C { w G D I \\{B* {ip, H'ijj)n){uj)\\u = a}, 
—av!{uj) = {B*{ip, H'ijj)-}{){uj), \u\-almost everywhere. 


Proof. Writing equation (26) as an integral yields 


(31) 

(32) 


+ {{B*{ip,H'ijj)n){uj),u'{oj))uj d|u|(a;) = 0. (33) 

For the integrand we obtain by the Cauchy-Bunyakovsky-Schwarz (CBS) inequality, using 


|u'(a;)|| = 1 and (27), 


Therefore (33) yields 


a + {{B*{ip,Hf;)-u){uj),u'{uj))u >a - \\{B*{ip, H'ijj)-u){uj)\\u > 0. 
a + {{B*{ip, Hfi))'u){oj),u {uj))u = 0 


(34) 


for |u|-almost all cu G D. For those oj the CBS inequality in (34) was sharp. This implies 
(32) and 

\\{B*{ip,Hfi:)'u){uj)\\u = a 

for lul-almost all a; G D. Since oj i—)• \\{B*{ip,H'if)'p){uj)\\u is continuous this implies 

d^l). □ 


For specific control operators B, equation (32) implies additional regularity for u'{oj). 
For example, in the case of the control operator from Example |3.1| we obtain the regularity 
u'{uj) G H^{<d,T). 

The relation (31) for the support of the optimal measure gives us the following corollary. 


Corollary 4.7. Let u be a local minimizer of (23). Then supp|tt| is compact. 

Proof. Since B*{(p, Hf;) G Co{Ll]U) we know that there is a compact set K C Ll such that 
B*{ip,H'4)) < a/2 for all uj ^ K. Using (31) this implies supplul C K. Therefore supplul 
is compact as a closed subset of the compact set K. □ 


This corollary is of significant physical interest. It says that although the frequency 
domain D might be unbounded, optimal solutions will always have bounded support. 
Controls from experiments, of course, always have this property, because arbitrarily fast 
oscillations are not realizable. But it is remarkable that such oscillations in our theoretical 
controls can be rules out rigorously. 

For most of the specific control operators considered in this work, a signihcantly stronger 
result holds. 
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Theorem 4.8 (Frequency sparsity). Let u he a loeal minimizer of (23) and let U and 


B be given as in one of the Examples 3.L 3.2. 3.3. Then supp|rt| is finite. 


Proof. The results follows from (31) by an analyticity argument. In the setting of the 


Examples 3.1-3.3 for every / G L'?(0, T), the map tv i—)• {B*f){uj) is real analytic. This is 


standard for the Fourier and Gabor transformation, and follows from analytic dependence 
on w of the right hand side in (18) for the dual two-scale operator. If 7^ M then B* f can 
be extended in the natural way to an analytic function on M. In all cases, {B* f)(u}) —?• 0 
if w —)■ ±00. 

For / = {(f, H'4 >)-h we obtain the analyticity of a; 1—)• || {B*{(p, i^Vj)'^)(a;)||^, which implies 
that the value is attained either everywhere in M or only at discrete points. Since the 
function vanishes at infinity and the support is compact by Corollary |4.7| the support 
condition (31) implies the finiteness of supp|tt|. □ 


This result implies that the optimal control is in fact a finite sum of Dirac measures 


in frequency with values in U. Example 3.4 cannot be treated with the same technique. 
Although the map u) 1 — {B* f ){oj) is complex analytic for / G L'^(0,T) if LI is identified 
with a subset of C, the analyticity does not carry over to to the absolute values. 


5 Numerical results 


In this section we apply the framework for sparse time-frequency control to different quan¬ 
tum systems. First we will describe our numerical approach. This includes a short dis¬ 
cussion of the discretization and the regularization of the optimal control problem. Then 
we will present two examples. The first example is the finite-dimensional system from 
Example 2.2 It serves to illustrate basic effects of sparse control of quantum systems. 
The second example is the two-level Schrodinger system from Example 2.3 The focus in 


this more challenging example will be the effect of different control spaces on the resulting 
optimal controls. 


5.1 Numerical approach 

Our numerical approach relies on the following three steps. Eirst, we discretize the mea¬ 
sure space by a finite sum of Dirac measures with values in a finite-dimensional Hilbert 
space. Then, we Huber-regularize the corresponding finite-dimensional nonsmooth prob¬ 
lem. Einally, we solve the resulting smooth optimization problem with a quasi-Newton 
method. 

The first step depends on D, of course. In our examples we always have D C 
for k G {1,2} and LI is an interval or the product of two intervals. We fix a uniform 
(tensor) grid and choose measures supported at those points as our ansatz space. Those 
measures can always be written as finite sums of Dirac measures. In this discrete setting 
the measure norm reduces to an norm for the coefficients from U multiplying the Dirac 
measures. To obtain a discrete problem we also need to discretize the Hilbert space U and 
the quantum system. In our examples we have lA = Hg(0,r;C) ovlA = L^(0, r;C) where 
we use piecewise linear finite elements on a uniform grid with the appropriate discrete 
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norms, or U = C were we do not need to discretize. The discrete Hilbert space is denoted 
by The control operator maps discrete controls to piecewise linear functions. The 
discretization of the quantum system depends on the system at hand. We approximate 
the time evolution of the discretized quantum system by a generalized Suzuki-Trotter 
method, see m- Together we obtain a finite-dimensional optimization problem that is 
non-smooth and non-convex. 

To deal with the nondifferentiability of the norm at the origin we Huber-regularize this 
non-smooth problem. We replace the norm of U in the norm by the function 

/i® : W ^ M given by 

, ^ / M\u - f, if \\z\\u > 0 , 

ImII^IIw> if||2||w<0, 

for some regularization parameter 9. The function h has the following two important 
properties: it is differentiable, and the derivatives of and the derivatives of the norm of 
U have the same behavior outside of a small neighborhood of zero. The first property makes 
the cost functional differentiable. The second property preserves the sparsity of solutions, 
in the sense that frequency regions with zero control amplitude turn into regions with 
control amplitude below the explicit threshold 9. More precisely we obtain the following 
theorem which holds e.g. in the case Q = [a;_,a;+] H /iZ for some mesh size h > 0. 


Theorem 5.1 (Characterization of frequency support after frequency discretization and 
Huber regularization). Let Q be a diserete set of admissable frequencies and choose any 
Huber regularization parameter 0 > 0. Consider the optimal control problem (12) with 
the norm replaced by the Huber regularized l^iU) norm Yluj&n h^{u{<^))- Then 

optimal controls u satisfy 


{w G n I \\u{uj)\\u >6'} = {a;Gfl| \\{B* {ip, H'ijj)){uj)\\u = a}. (35) 


Proof. In this case the cost is differentiable and the control is a function rather than a 
measure with respect to frequency, and the first order optimality condition 0 = j'{u){5u) 
acquires the simple form 

0 = {B* {ip, H'ip)){u!) +aVh^{u{uj)) for all io £ Q. 

Since Vh^{z) equals zj^z^u for || 2 ;|| > 9 and zl9 for \\z\\u < 9, the assertion follows. □ 

In the numerical examples below, only a few frequencies were contained in the superlevel 
set on the left which numerically replaces the frequency support of u. See e.g. Figure 

We solved the resulting smooth problem with a quasi-Newton method. Since the di¬ 
mension of the control space can become quite large with our approach we chose the 
memory efficient L-BFGS method, see [2^. For the numerical realization gradients for the 
discretized problems were used. The optimization method was terminated as soon as the 
norm of the gradient relative to the largest gradient was below a tolerance of at 
least 10“®. 

The cost parameter is always chosen such that one achieves at least 95% of the control 
objective, i.e. ^{'i/j{T),0'ip{T))'p < 0.025. An automated choice of a would be helpful to 
obtain useful cost parameters for a variety of problems. 
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Figure 2: Optimal control when the cost is chosen as a measure norm with respect to 
frequency and the Hq norm with respect to time (Example |2.2[ ). (a) Control held 
as a function of time, (c) Time-frequency representation u{uj, t) (color indicates 
absolute value), (b) The contributions due to the two active frequencies of the 
optimal held. 


The result of this nonlinear optimization problem also depends on the initial guess for 
the control. For small a, the initial guess for the control uses a hxed element in together 
with a random complex phase for the different frequencies. For larger a, where snch a 
generic initial guess leads to convergence of the method to suboptimal critical points near 
the origin, we nse optimal solutions for smaller a as initial guess. 


5.2 Three level system 

As our hrst example we chose a typical hnite-dimensional projection of an atom in a laser 


held, see Example 2.2 We use the matrices 



(-2 

0 

0\ 

1 j 

(0 

0 

1\ 

II 

0 

-1 

0 

, Hi =\ 

0 

0 

1 


^ 0 

0 

2 

1 1 


1 

0/ 


This corresponds for instance to a Is, 2s, and 3p state. Note that the coupling matrix 
elements Q between the Is and 2s state vanishes, i.e. the transition Is —?• 2s is “forbidden” 
since these states have equal parity. On the other hand the transitions Is —)• 3p and 3p 
—>■ 2s are allowed. The control objective is to reach the third eigenstate starting from 
the ground state. The initial condition and the observation operator are then given by 
'ipo = (1,0,0) and O = diag(l, 1,0). We use a frequency band 0 = [2,5] discretized with 
100 grid points. The expected transition frequencies (Ui = 3 and a;2 = 4 are contained 
in n. We chose a time horizon of T = 100 and a time grid with 4096 points. The time 
horizon was chosen large enough to allow for sufficiently many periods with the transition 

l«IU(n;Hi(o,T;C)) with the control 


3.1 


i.e. 


freqnencies. We chose the cost from Example 
operators given by @-([17]), and a cost parameter of a = 0.1. The control objective was 
achieved to more than 99.999%. 

The great advantage of the measure space control is that we can decompose the field into 
simple components. Figure]^ shows the optimal control. We see that only two frequencies 
have a visible contribution. They correspond to the two Bohr frequencies wi and a;2- The 


18 




































frequency frequency 

(a) (b) 


Figure 3: Detailed numerical illustration of frequency sparsity of the optimal control from 
Figure]^ (b) By Theorem 5.1 the optimal control should vanish off the frequen¬ 
cies where the norm \\{B*{ip, H’4})){u})\\u (dots in (b)) reaches the cost parameter 
a (solid line in (b)). (a) The numerical optimal control (dots in (a)), instead of 
vanishing, drops by three orders of magnitude below the threshold given by the 
numerical regularization parameter 6 (dashed line in (a)), precisely as theoreti¬ 
cally predicted by equation (35). 


time profiles are smoothly switched on and off and remain active during the whole time. 
This is consistent with the choice U = i?g(0, T;C), which promotes smoothness and non¬ 
locality in time. The field for the first transition reaches its maximum before the field for 
the second transition. This reflects the intuitive idea that we have to induce the transition 
between levels one and three before that between levels three and two. 

In fact, unlike most previous control-theory-based forcing fields (for an exception see |23| 
where a weighted regularization is used), the field obtained here bears considerable 
resemblance to the simple and intuitive few-parameter pulses which have been used by 
laser physicists for a long time. Compare, in particular, the two pulses in Figure 9 of [1], 
whose achievement of the control objective was beautifully demonstrated experimentally. 

Figure gives a more detailed numerical illustration of the frequency sparsity of the 
optimal control. Due to the effect of Huber regularization with 0 > 0, the numerical 
analog of the support of u is the superlevel set {w | ||u(a;)||w > & } (see equation (35)), 
which is seen to consist of just two frequencies. 


5.3 Schrodinger dynamics on two potential energy surfaces 

In this second example we consider a Schrodinger system on two potential energy surfaces 


as arising in the laser control of chemical reactions, see Example 2.3 The spectral gap 


between the two potential energy surfaces varies depending on the position of the nuclear 
wave function and therefore a much larger variety of frequencies is potentially useful for 
successful control. We also expect an additional time structure in the controls due to 
the movement of the densities in space. Therefore it is much more challenging to obtain 
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simple controls for this example compared to the previous one. We focus on a comparison 
between controls generated for different choices of the space lA of time profiles and the 
control operator B. 

For simplicity we limit ourselves to one reaction coordinate, i.e. space dimension d = 1. 
The control objective is to reach the potential well on the right starting from the potential 
well on the left, see Figure The initial state ipo is a Gaussian located in the lower well. 
The observation operator O is the projection on the complement of functions with support 
on the lower surface on the right of the potential barrier. Instead of the physical coupling 
given by (10) we use a coupling operator given by the multiplication with the reaction 
coordinate on the diagonal and the identity on the off-diagonal. 

The energy differences between the two potential energy surfaces measured at the local 
minima of the lower surface are approximately 0.074 and 0.048. Therefore we expect 
optimal controls to contain the two frequencies oji k. 0.074 and W 2 ~ 0.048. We chose a 
time horizon of T = 3000 and a time grid with 2048 points. The time horizon was chosen 
large enough to allow for sufficiently many periods with the Bohr frequencies wi and a; 2 , 
and for sufficient movement of the wave function in space. For the discretization in space 
we use a simple finite difference scheme on the domain [—4,4] with 256 grid points. We 
validated that this discretization is sufficient for the range of parameters relevant to this 
application. 

We compare the resulting optimal fields for different choices of the frequency domain 
Q, the Hilbert space lA of admissible time profiles, and the control operator B. We also 
compare them to the optimal field for the classical Hilbert space cost functional with 
L^(0,T;M) norm. In particular we choose the following setups. 


Two-scale synthesis, smooth functions of time 

HQ(0,r;C), B as in Example 3.1 


n = [1/30,1/10], u = 
The frequency band contains the expected 
transition frequencies oJi and a; 2 - It is discretized with 100 grid points. We dis¬ 
cretize lA with linear finite elements. The time grid corresponds to the grid of the 
time stepping. This results in 100 • 2 • 2048 = 409600 real degrees of freedom. 


• Gabor synthesis, functions of time. H = [1/30,1/10], U = L^(0,r;C) 
and B as in Example |3.2[ with the Gaussian kernel k suitably adapted to generate 
homogeneous Dirichlet boundary conditions. We discretize 12 and lA as before. Eor 
the evaluation of B we explicitly construct the matrix K corresponding to the kernel 
k. 


• Fourier synthesis, constant functions of time. Q = [1/30,1/10], lA = C, B the 
Eourier synthesis operator, see Example |3.3[ The frequency band 12 is discretized 
with 100 grid points. This results in 2 • 100 = 200 real degrees of freedom. 

• Time-frequency Gabor synthesis. 12 = [1/30,1/10] x [0,T], lA = C, B as in 
Example |3.4[ The time-frequency cylinder 12 is discretized by a tensor grid. In 
frequency direction we use a regular grid with 100 grid points. In time direction we 
use a grid of 14 points. This results in 100 • 14 • 2 = 2800 real degrees of freedom. 

• Standard cost. This means we directly minimize over v G L^(0,r;M) the 
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funtional J(V',f) = \{'il^{T),0'ip{T))) + a\\v\\\^ 2 - We use linear finite elements on the 
time grid of the time stepping method. This results in 2048 real degrees of freedom. 


For the simulations we chose a such that probability to end up in the desired subspace 
is near the value of 95%. 

Figurej^shows the optimal controls. We see that their structure depends on the different 
cost functionals, and that the measure cost (13) for the two-scale and Gabor time-frequency 
representation leads to an extremely small number of active frequencies. In the following 
we will discuss the most signihcant differences and similarities. 

Turning to the second column of Figure]^ we focus our attention to the regions around 
the two Bohr frequencies coi and UJ 2 - The regions correspond to the transitions up from the 
hrst well and down into the second well, respectively. The expected and desired behavior 
of the frequencies is clearly obtained in the hrst three choices of the cost functional. 
This structure becomes less pronounced for the time-frequency Gabor approach, where 
the stronger localization of the pulses in time leads to broader regions in the frequency 
representation, and for the control, which has the frequency prohle that is most difficult 
to interpret. 

The third column of Figure illustrates the sparse time-frequency structure of the 
optimal controls. To compare with the standard control v, we also computed a 
sparse representation u of the latter, by a posteriori minimization of \\Bu — '^ 11 ^ 2 ( 0 '^) + 


®ll'^llA4(o //i(o,'r;C)) with the two-scale synthesis operator B and a = 10“^. 

Figurej^shows that time-frequency representation of the control contains many more 

active frequencies compared to the other controls. 

It is interesting that the space M. (fl; Hq (0, T; C)) leads to more active frequencies com¬ 
pared to the approach with space A4(0; L^(0, T; C)) with the Gabor synthesis operator. 
The Fourier and the time-frequency Gabor approaches lead to larger but still reasonably 
small number of active frequencies. The choice of a specihc cost functional may be guided 
by the concrete application under consideration. 


6 Conclusion 

In summary, measure valued costs imposed on time-frequency representations of the elec¬ 
tric held as introduced in this paper systematically produce frequency sparse controls, in 
contrast to standard costs. The resulting controls resemble physically intuitive controls 
designed by experimentalists. We hope that the measure-space, time-frequency approach 
will reduce the current gap between numerical controls on the one side and experimental 
implementation and physical intuition on the other. The hexibility of the approach can 
be further exploited to construct controls suited for concrete experimental setups. 

Acknowledgements. The work of F. Henneke was funded by DFG through the Inter¬ 
national Research Training Group IGDK 1754. 
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Figure 4: Optimal controls for Schrodinger dynamics on two potential energy surfaces. 

Rows; Different choices of cost functionals and control operators. Columns: 
Time, frequency, and time-frequency representation of the optimal controls (i.e. 
{Bu){t), Bu{t) and u(a;,t)). In the rightmost column, the absolute values of the 
optimal measures are plotted in the time-frequency plane. 


22 























































































































References 


[1] A. Auger, A. Ben, Haj Yedder, E. Cances, C. Le Bris, C. M. Dion, A. Keller, and 
O. Atabek. Optimal laser control of molecular systems: Methodology and results. 
Math. Models Methods Appl. Sei, pages 1281-1315. 

[2] G. G. Balint-Kurti, S. Zou, and A. Brown. Optimal Control Theory for Manipulating 
Molecular Processes., pages 43-94. John Wiley &: Sons, Inc., 2008. 

[3] J. M. Ball, J. E. Marsden, and M. Slemrod. Gontrollability for distributed bilinear 
systems. SIAM J. Control Optim., 20(4);575-597, 1982. 

[4] H. Bergmann, H. Theuer, and B. W. Shore. Goherent population transfer among 
quantum states of atoms and molecules. Rev. Mod. Phys., 70(3): 1003-1025, 1998. 

[5] U. Boscain, G. Gharlot, J.-P. Gauthier, Stephane Guerin, and Hans-Rudolf Jauslin. 
Optimal control in laser-induced population transfer for two- and three-level quantum 
systems. J. Math. Phys.., 43(5):2107-2132, 2002. 

[6] E. Gasas, R. Herzog, and G. Wachsmuth. Optimality conditions and error analysis 
of semilinear elliptic control problems with cost functional. SIAM J. Optim., 
22(3):795-820, 2012. 

[7] T. Ghambrion. Periodic excitations of bilinear quantum systems. Automatica J. 
IFAC, 48(9):2040-2046, 2012. 

[8] F. H. Glarke. Optimization and nonsmooth analysis, volume 5 of Classics in Applied 
Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, 
PA, second edition, 1990. 

[9] Ghristian Glason and Karl Kunisch. A duality-based approach to elliptic control prob¬ 
lems in non-reflexive Banach spaces. ESAIM Control Optim. Calc. Var., 17(1):243- 
266, 2011. 

[10] D. D’Alessandro. Introduction to quantum control and dynamics. Chapman 
& Hall/CRC Applied Mathematics and Nonlinear Science Series. Chapman &: 
Hall/CRC, Boca Raton, EL, 2008. 

[11] Steffen J Glaser, Ugo Boscain, Tommaso Calarco, Christiane P Koch, Walter 
Kockenberger, Ronnie Kosloff, Ilya Kuprov, Burkhard Luy, Sophie Schirmer, Thomas 
Schulte-Herbriiggen, et al. Training schrodingers cat: Quantum optimal control. The 
European Physical Journal D, 69(12):l-24, 2015. 

[12] M. Hellgren, E. Rasanen, and E. K. U. Gross. Optimal control of strong-held ion¬ 
ization with time-dependent density-functional theory. Phys. Rev. A, 88:013414, Jul 
2013. 

[13] F. Henneke and M. Liebmann. A generalized Suzuki-Trotter type method in optimal 
control of coupled Schrodinger equations, pre-print, 2015. 


23 



[14] R. Herzog, G. Stadler, and G. Wachsmuth. Directional sparsity in optimal control of 
partial differential equations. SIAM J. Control Optim., 50(2):943-963, 2012. 

[15] M. Hintermiiller, D. Maralirens, P. A. Markowicli, and G. Sparber. Optimal bilinear 
control of Gross-Pitaevskii equations. SIAM J. Control Optim., 51(3):2509-2543, 
2013. 

[16] P. von den Hoff, S. Thallmair, M. Kowalewski, R. Siemering, and R. de Vivie-Riedle. 
Optimal control theory - closing the gap between theory and experiment. Phys. Chem. 
Chem. Phys., 14:14460-14485, 2012. 

[17] K. Ito and K. Kunisch. Optimal bilinear control of an abstract Schrodinger equation. 
SIAM Journal on Control and Optimization, 46(l);274-287, 2007. 

[18] K. Kormann, S. Holmgren, and H.O. Karlsson. A Fourier-coefficient based solution 
of an optimal control problem in quantum chemistry. Journal of Optimization Theory 
and Applications, 147:491-506, 2010. 

[19] K. Kunisch, K. Pieper, and B. Vexler. Measure valued directional sparsity for 
parabolic optimal control problems. SIAM J. Control Optim., 52(5):3078-3108, 2014. 

[20] S. Lang. Real and functional analysis, volume 142 of Craduate Texts in Mathematics. 
Springer-Verlag, New York, third edition, 1993. 

[21] M. Lapert, D. Sugny, R. Tehini, and G. Turinici. Monotonically convergent optimal 
control theory of quantum systems with spectral constraints on the control field. 
Physical Review A: Atomic, Molecular and Optical Physics, 79:063411, 2009. 

[22] X. J. Li and J. M. Yong. Optimal control theory for infinite-dimensional systems. 
Systems &: Gontrol: Foundations &: Applications. Birkhauser Boston Inc., Boston, 
MA, 1995. 

[23] Yvon Maday, Julien Salomon, and Gabriel Turinici. Monotonic parareal control for 
quantum systems. SIAM J. Numer. Anal, 45(6):2468-2482 (electronic), 2007. 

[24] L. Meziani. On the dual space (70(5, A). Acta Math. Univ. Comenian. (N.S.), 
78(1):153-160, 2009. 

[25] J. Nocedal and S. J. Wright. Numerical optimization. Springer Series in Operations 
Research and Financial Engineering. Springer, New York, second edition, 2006. 

[26] A. P. Peirce, M. A. Dahleh, and H. Rabitz. Optimal control of quantum-mechanical 
systems: existence, numerical approximation, and applications. Phys. Rev. A (3), 
37(12):4950-4964, 1988. 

[27] 1. 1. Rabi. Space quantization in a gyrating magnetic field. Phys. Rev., 51:652-654, 
Apr 1937. 


24 



[28] Michael Reed and Barry Simon. Methods of modern mathematical physics. II. Fourier 
analysis, self-adjointness. Academic Press [Harcourt Brace Jovanovich, Publishers], 
New York-London, 1975. 

[29] Q. Ren, G. G. Balint-Kurti, F. R. Manby, M. Artamonov, T.-S. Ho, and H. Rabitz. 
Quantum control of molecular vibrational and rotational excitations in a homonuclear 
diatomic molecule: A full three-dimensional treatment with polarization forces. The 
Journal of Chemical Physics, 124(1):014111, 2006. 

[30] S. Ruetzel, C. Stolzenberger, F. Dimler, D. J. Tannor, and T. Brixner. Adaptive 
coherent control using the von Neumann basis. Phys. Chem. Chem. Phys., 13:8627- 
8636, 2011. 

[31] J. Scheuer, X. Kong, R. S. Said, J. Chen, A. Kurz, L. Marseglia, J. Du, P. R. Hemmer, 
S. Montangero, T. Calarco, B. Naydenov, and F. Jelezko. Precise qubit control beyond 
the rotating wave approximation. New Journal of Physics, 16(9):093022, 2014. 

[32] S. Sharma, H. Singh, and G. G. Balint-Kurti. Genetic algorithm optimization of laser 
pulses for molecular quantum state excitation. The Journal of Chemical Physics, 
132(6):064108, 2010. 

[33] B. W. Shore. The theory of coherent atomic excitation. Wiley, New York, NY, 1990. 

[34] G. Stadler. Elliptic optimal control problems with L^-control cost and applications 
for the placement of control devices. Comput. Optim. AppL, 44(2):159-181, 2009. 

[35] S. Teufel. Adiabatic perturbation theory in quantum dynamics, volume 1821 of Lecture 
Notes in Mathematics. Springer-Verlag, Berlin, 2003. 

[36] G. Turinici, C. Le Bris, and H. Rabitz. Efficient algorithms for the laboratory dis¬ 
covery of optimal quantum controls. Phys. Rev. E, 70:016704, Jul 2004. 

[37] G. von Winckel and A. Borzi. Computational techniques for a quantum control 
problem with H^-cost. Inverse Problems, 24(3):034007, 2008. 

[38] G. von Winckel, A. Borzi, and S. Volkwein. A globalized Newton method for the 
accurate solution of a dipole quantum control problem. SIAM J. Sci. Comput., 
31(6):4176-4203, 2009/10. 

[39] G. Vossen and H. Maurer. On L^-minimization in optimal control and applications 
to robotics. Optimal Control Appl. Methods, 27(6):301-321, 2006. 


25 



